CHIRAL SYMMETRY BREAKING BY INSTANTONS ft 



q\ ! Dmitri Diakonov 



x 



Petersburg Nuclear Physics Institute, Gatchina, St. Petersburg 188350, Russia 



m 

> 

! 1 Foreword 

m 

<N . 

O Instantons are certain configurations of the Yang-Mills potentials A a (x) satisfying 

the equations of motion D a ^F^ u = in euclidean space, i.e. in imaginary time. The 
^ solution has been found by Belavin, Polyakov, Schwartz and Tiupkin in 1975 ffl; the 

name "instanton" has been suggested by 't Hooft in 1976 ||, who also made a major 
D ! contribution to the investigation of the instantons properties. 

In QCD instantons are the best studied non-perturbative effects, leading to the 
formation of the gluon condensate and of the so-called topological susceptibility 
needed to cure the U(l) paradox 0, |J. The QCD instanton vacuum has been studied 
starting from the pioneering works in the end of the seventies |5|, |6|; a quantitative 
treatment of the instanton ensemble, based on the Feynman variational principle, has 
been developed in ref. 0. The most striking success of the QCD instanton vacuum 
is its capacity to provide a beautiful mechanism of the spontaneous chiral symmetry 



breaking || || 10]. Moreover, the instanton vacuum leads to a very reasonable effec- 
tive chiral lagrangian at low energies, including the Wess-Zumino term, etc., which, 
in its turn gives a nice description of nucleons as chiral quark solitons [|ll|] . 

It should be stressed that literally speaking instantons do not lead to confinement, 
though they induce a growing potential for heavy quarks at intermediate separations 
[fl^l ; asymptotically it flattens out [|, However, it has been realized a decade 
ago ^ [nj], that it is chiral symmetry breaking and not confinement that deter- 
mines the basic characteristics of light hadrons (one would probably need an explicit 
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confinement to get the properties of highly excited hadrons). Therefore, since the 
instanton vacuum describes well the physics of the chiral symmetry breaking, one 
would expect that instantons do explain the properties of light hadrons, both mesons 
and baryons. Indeed, a detailed numerical study of dozens of correlation functions in 
the instanton medium undertaken by Shuryak, Verbaarschot and Schafer [15| (earlier 
certain correlation functions were computed analytically in refs. || ||) demonstrated 
an impressing agreement with the phenomenology |16| and, more recently, with di- 
rect lattice measurements [jDJ. As to baryons, the instanton-motivated chiral quark 
soliton model [jlTJ also leads to a very reasonable description of dozens of baryon 



characteristics (for a review see ref . [jig] ) . 

More recently the instanton vacuum was studied in direct lattice experiments by 



the so-called cooling procedure [19], ^y, |2l], g^, [17], Pq] , see also the proceedings of the 
Lattice-94 meeting for a more complete list of references. It was demonstrated that 
instantons and antiinstantons (J's and J's for short) are the only non-perturbative 
gluon configurations surviving after a sufficient smearing of the quantum gluon fluc- 



tuations. The measured properties of the II medium appeared |21|, [17], ^3] to be close 
to that computed from the variational principle [0] and to what had been suggested 
by Shuryak in the beginning of the 80's || from phenomenological considerations. 

Cooling down the quantum fluctuations above instantons kills both the one-gluon 
exchange and the linear confining potential (a small residual string tension observed 



in the cooled vacuum |2T], [L7| is probably due to the instanton-induced rising potential 
at intermediate distances [0). Nevertheless, in the cooled vacuum where only J's 
and J's are left, the correlation functions of various mesonic and baryonic currents, as 
well as the density-density correlation functions appear to be quite similar to those of 
the true or "hot" vacuum JT7J. We consider it to be a remarkable confirmation of the 
ideology which has been put forward quite some time ago @, ||, || [14J : instantons are 
responsible for the basic properties of the QCD vacuum and of the low-mass hadrons, 
while the confinement, fundamental as it is, does not readily manifest itself in those 
properties. 

These lectures are aimed as an introduction to understanding the above works. 
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2 Periodicity of the Yang— Mills potential energy 



I start by explaining the physical meaning of instantons as classical tunneling tra- 
jectories in imaginary (euclidean) time. To that end I shall temporarily work in the 
Aq = gauge, called Weyl or Hamiltonian gauge, and forget about fermions for a 
while. The remaining pure Yang-Mills or "pure glue" theory is nonetheless non- 
trivial, since gluons are self-interacting. For simplicity I start from the SU(2) gauge 
group. 

The spatial YM potentials Af(pc,t) can be considered as an infinite set of the 
coordinates of the system, where i — 1, 2, 3, a = 1, 2, 3 and x are "labels" denoting 
various coordinates. The YM action is 

s = i? I d ' x = S dt (w I d3x E2 - h f** B2 ) ( "> 

where E is the electric field stregth, 

E? (x, t) = A^x, t) (2.2) 
(the dot stands for the time derivative), and B is the magnetic field stregth, 

Bf (x, t) = \e ijk (djAt - d k A) + e^A'At) . (2.3) 

Apparently, the first tirm in eq. ( |2.1| ) is the kinetic energy of the system of co- 
ordinates {Af(pc, t)} while the second term is minus the potential energy being just 



the magnetic energy of the field. The simple and transparent form of eq. (2.2) is the 
advantage of the Weyl gauge. 

Let us introduce an important quantity called the Pontryagin index or the four- 
dimensional topological charge of the YM fields: 

If 1 

Qt = J d 4 x Fj^F^,, F* v = -e^apF^p. (2.4) 

The integrand in eq. ( |2.4|) is a full derivative of a four- vector K^. 

^ F ^v = K » = i^a Wr (Kd,A; + \e abc AlA h p A^ . (2.5) 

Therefore, assuming the fields A^ are decreasing rapidly enough at spatial infinity, 
one can rewrite the 4-dimensional topological charge (|2.4|) as 
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Qt = J d 4 x{d K - diKi) = Jdt-J d 3 xK . (2.6) 
Introducing the Chern-Simons number 

N cs = Jd 3 xK = ^J rf 3 x e ijk {A^Al + l -e abc A«A)Al) (2.7) 

we see from eq. ( |2.6|) that Qt can be rewritten as the difference of the Chern-Simons 
numbers characterizing the fields at t — ±00: 

Qt = N cs (+oo) - N cs (-oo). (2.8) 

The Chern-Simons characteristics of the field has an important property that it 
can change by integers under large gauge transformations. Indeed, under a general 
time-independent gauge transformation, 

Ai-ttfAiU + irfdiU, Ai = A""7T> (2-9) 
the Chern-Simons number transforms as follows: 

N CS -> N cs + N w (2.10) 



where N\y is the winding number of the gauge transformation ( |2.9|) : 

Nw= 2^I ^ X €ijk [( f/t ^ f/ )( f/t ^ f/ )( f/t ^ f/ )] • ( 2 - n ) 
The SU(2) unitary matrix U of the gauge transformation ( |2.9| ) can be viewed as 
a mapping from the 3-dimensional space onto the 3-dimensional sphere of parameters 
S3. If at spatial infinity we wish to have the same matrix U independently of the 
way we approach the infinity (and this is what is usually assumed), then the spatial 
infinity is in fact one point, so the mapping is topologically equivalent to that from 
S3 to S3. This mapping is known to be non-trivial, meaning that mappings with 
different winding numbers are irreducible by smooth transformations to one another. 
The winding number of the gauge transformation is, analytically, given by eq. Q2.ll 



As it is common for topological characteristics, the integrand in ( |2.11| ) is in fact 
a full derivative. For example, if we take the matrix U(x) in a "hedgehog" form, 
U = exp[i(rr/r ■ P(r)], eq. (|2.11| ) can be rewritten as 
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N w = - I dr sm 2 P— = - 

dr Ti 



7T 



sin 2P 



= integer (2-12) 

o 



since P(r) both at zero and at infinity needs to be multiples of it if we wish U (r) to 
be unambigiously defined in the origin and at the infinity 
Let us return now to the potential energy of the YM fields, 

V=^jd**{B«f. (2.13) 

One can imagine plotting the potential energy surfaces over the Hilbert space of 
the coordinates A"(x). It will be some complicated mountain country, like around 
Varenna. If the field happens to be a pure gauge, Ai = iWdiU, the potential energy 
at such points of the Hilbert space is naturally zero. Imagine that we move along 
the "generalized coordinate" being the Chern-Simons number fl2.7|) , fixing all other 
coordinates whatever they are. Let us take some point Af(x) with the potential 
energy V. If we move to another point which is a gauge transformation of Af(x) with 
a winding number N\y, its potential energy will be exactly the same as it is strictly 
gauge invariant. However the Chern-Simons "coordinate" of the new point will be 
shifted by an integer Nw from the original one. We arrive to the conclusion first 
pointed out by Faddeev and Jackiw and Rebbi [^5] in 1976, that the potential 
energy of the YM fields is periodic in the particular coordinate called the Chern- 
Simons number. 

One may wish to plot this periodic dependence of the YM potential energy on 
Ncs ■ Putting such a problem in a situation where the potential energy depends also 
on an infinite number of other " coordinates" , we imply that one is to find the minimal 
energy path, say, from N C s = to N C s = 1 , for a given value of N C s ■ However, 
the pure YM theory is scale-invariant at the classical level, so to solve the problem 
one has to fix the spatial size of the Ai fields somehow. The situation is different in 
the electro- weak theory where the scale invariance is explicitly broken at the classical 
level by the non-zero Higgs vacuum expectation value. Therefore the problem of 
finding the minimal energy pass in the mountain country of the YM Hilbert space 
is well defined. It was solved by Akiba, Kikuchi and Yanagida in 1988 [p6|1 , and the 
reader may satisfy his or her curiosity and have a view of the periodic dependence 
of the potential energy on the Chern-Simons number in that work (see also ref, 
where this dependence is generalized to non-zero matter density). 
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The static configuration of the A{ fields corresponding to the saddle point in the 
minimal-energy path and having exactly Ncs = 1/2 is called sphaleron. It was found 
many years ago by Dashen, Hasslacher and Neveu |28[ and rediscovered in the context 



of the electro- weak theory by Manton and Klinkhammer |29[. The sphaleron is the 
"coordinate" of the top of the potential barrier separating two zero-potential points, 
Ncs — and Ncs = 1 • In the electro- weak theory the height of the barrier is of 
the order of my/fa {my/ is the mass of the W boson and a is the gauge coupling 
constant). In an unbroken YM theory like QCD the classical energy barrier between 
topologically distinct vacua can be made infinitely small due to the scale invariance. 
However that does not mean that the barriers are easily penetrable: we shall calculate 
the transition amplitudes from one mimimum to another in the next section. 



3 Instanton configurations 

In perturbation theory one deals with zero-point quantum-mechanical fluctuations of 
the YM fields near one of the minima, say, at Ncs — • The non-linearity of the YM 
theory is taken into account as a perturbation, and results in series in g 2 where g is 
the gauge coupling. In that approach one is apparently missing a possibility for the 
system to tunnel to another minimum, say, at Ncs = 1 . The tunneling is a typical 
non-perturbative effect in the coupling constant, and instantons have direct relation 
to the tunneling. 

The tunneling amplitude can be estimated as exp(-S), where S is the action along 
the classical trajectory in imaginary time, leading from the minimum at Ncs — at 
t = — oo to that at Ncs = 1 at £ = +oo f3~0| . According to eq. ( |2.8| ) the 4-dimensional 
topological charge of such trajectory is Qt = 1. To find the best tunneling trajectory 



having the largest amplitude one has thus to minimize the YM action (|2.1|) provided 
the topological charge ( |2.4|) is fixed to be unity. This can be done using the following 
trick |]. Consider the inequality 



o < / d 4 x (f; v - i> 

J d 4 x (2F 2 - 2FF) = 8g 2 S - 64vr 2 Q T , (3.1) 



hence the action is restricted from below: 
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S> S 4Qt = K- (3-2) 

g 2 g 2 

Therefore, the minimal action for a trajectory with a unity topological charge is 
equal to 87r 2 /g 2 , which is achieved if the trajectory satisfies the self-duality equation: 

F^ u = F{lv (3-3) 

Notice that any solution of eq. ( p^3|) is simultaneously a solution of the general 
YM equation of motion D a ^F^ u = 0: that is because the "second pair" of the Maxwell 
equations, D^F^ U = 0, is satisfied identically 

To solve eq. (|3~3j ) let us recall a few facts about the Lorentz group SO(3, 1). Since 
we a talking about the tunneling fields which can only develop in imaginary time, it 
means that we have to consider the fields in euclidean space-time, so that the Lorentz 
group is just 5*0(4) = SU(2) x SU(2). The gauge potentials belong to the (|, |) 
representation of the SU(2) x SU{2) group, while the field strength F^ v belongs to 
the (1,1) representation. In other words it means that one linear combination of F^ v 
transforms as a vector of the left SU(2), and another combination transforms as a 
vector of the right SU(2). These combinations are 

Ft = VU F ^ + *U f r = VU F ^ ~ 3-0. ( 3 - 4 ) 
where 77,77 are the so-called 't Hooft symbols described in ref.[Q, see also below. We 
see therefore that a self-dual field strength is a vector of the left SU(2) while its right 
part is zero. Keeping that experience in mind we look for the solution of the self-dual 
equation in the form 

^ = Vva X v ""■ (3-5) 

Using the formulae for the 77 symbols from ref.[H one can easily check that the 
YM action can be then rewritten as 



1 fd§\ 2 1 



+ -($ 2 - l) 2 



in i 7 . (3.6) 



2\dr 

This can be recognized as the action of the double-well potential whose minima 
lie at $ = ±1, and r plays the role of "time", p is an arbitrary scale. The trajectory 
which tunnels from 1 at r = — 00 to— 1 at r = +00 is 
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$ = -tanh(£), (3.7) 
and its action ( |3.6[ ) is S = 8ir 2 /g 2 , as needed. Substituting the solution ( |3.7[ ) into 



31) we get 



w = jy., - ( 3 - 8 ) 



^ v ' x 2 (x 2 + p 2 ) 
The correspondent field strength is 



F a = ^ (V - 2^ XaXv - 2n p X » X P\ F a F a = 192 p4 (3 9) 

^ (x 2 + p 2 ) 2 V ha laa x 2 lva x 2 /' » v ^ (x 2 + p 2 ) 4 ' y } 



and satisfies the self-duality condition (|3.3p . 

The anti-instanton corresponding to tunneling in the opposite direction, from 
N C s = 1 to N C s = , satisfies the antz-self-dual equation, with F — > —F; its 



concrete form is given by eqs. (|3.8| , |3.9|) with the replacement 77 — > 77 



Eqs.(|3~8|,[3.9|) describe the field of the instanton in the singular Lorentz gauge; the 



singularity of at x = is a gauge artifact: the gauge- invariant field strength 
squared is smooth at the origin. The formulae for instantons are more simple in the 
Lorentz gauge, and I shall use it further on. 

The instanton field, eq. (|3.8|) , depends on an arbitrary scale parameter p which we 
shall call the instanton size, while the action, being scale invariant, is independent of 
p. One can obviously shift the position of the instanton to an arbitrary 4-point - 
the action will not change either. Finally, one can rotate the instanton field in colour 
space by constant unitary matrices U. For the SU{2) gauge group this rotation is 
characterized by 3 parameters, e.g. by Euler angles. For a general SU(N C ) group the 
number of parameters is N 2 — 1 (the total number of the SU(N C ) generators) minus 
(N c — 2) 2 (the number of generators which do not effect the left upper 2x2 corner 
where the standard SU(2) instanton ( |3.8| ) is residing), that is 4iV c — 5. These degrees 
of freedom are called instanton orientation in colour space. All in all there are 

4 (centre) + 1 (size) + (4iV c — 5) (orientations) = 4N C (3.10) 

so-called collective coordinates desribing the field of the instanton, of which the action 
is independent. 
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It is convenient to indroduce 2x2 matrices 



a± = (±ilt,l), x ± =x,a±, (3.11) 

such that 

^r a V ^^a+a;-a+a-; 2ir a f% a = a~a+ - a~a+, (3.12) 

then the instanton field with arbitrary center z^, size p and colour orientation U in 
the SU (N c ) gauge group can be written as 



-ip 2 U\a-(x - z) + - (x- z) u ]Ul h 1 h 

or 



?n 2 (l ab fi fJ ' (r — A 

A" = ±£ ''^ bK , } * , O ab = Tr(UH a Ua b ), O ab O ac = 5 bc . (3.14) 
F (x — z) 2 [p 2 + (x — z) 2 \ 

Physically, one can think of instantons in two ways: on one hand it is a tunneling 
process occuring in imaginary time (this interpretation belongs to V.Gribov, 1976), 
on the other hand it is a localized pseudoparticle in the euclidean space (A.Polyakov, 



1977 [3l|) 



4 Gluon condensate 

The QCD perturbation theory implies that the fields A^(x) are performing quantum 
zero-point oscillations; in the lowest order these are just plane waves with arbitrary 
frequences. The aggregate energy of these zero-point oscillations, (B 2 + E 2 )/2, is 
divergent as the fourth power of the cutoff frequency, however for any state one has 
(F 2 ) = 2(B 2 — E 2 ) = 0, which is just a manifestation of the virial theorem for 
harmonic oscillators: the average potential energy is equal to that of the kinetic (I 
am temporarily in the Minkowski space). One can prove that this is also true in 
any order of the perturbation theory in the coupling constant, provided one does 
not violate the Lorentz symmetry and the renormalization properties of the theory. 
Meanwhile, we know from the QCD sum rules phenomenology that the QCD vacuum 
posseses what is called gluon condensate B: 
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3^2 (CO - (200 MeVf > 0. (4.1) 

Instantons suggest an immediate explanation of this basic property of QCD. In- 
deed, instanton is a tunneling process, it occurs in imaginary time; therefore in 
Minkowski space one has Ef = ±iBf (this is actually the duality eq. (|3.3|) ). There- 
fore, one gets a chance to explain the gluon condensate. In euclidean space the electric 
field is real as well as the magnetic one, and the gluon condensate is just the average 
action density. Let us make a quick estimate of its value. 

Let the total number of J's and J's in the 4-dimensional volume V be N. Assuming 
that the average separations of instantons are larger than their average sizes (to be 
justified below), we can estimate the total action of the ensemble as the sum of 
invidual actions (see eq. ( p.2[) ): 

{rf*)V = J d^xFl ~ N ■ 327T 2 , (4.2) 

hence the gluon condensate is directly related to the instanton density in the 4- 
dimensional euclidean space-time: 

1 N 1 

~ * T7 = -m- (4-3) 



32tt 2X ^ V R 4 ' 

In order to get the phenomeno logical value of the condensate one needs thus to have 
the average separation between pseudoparticles |J 

R ~ 1 = 1 fm. (4.4) 

200 MeV V ; 

There is another point of view on the gluon condensate which I describe briefly. 
In principle, all information about field theory is contained in the partition function 
being the functional integral over the fields. In the euclidean formulation it is 



Z = j DA.exp (-± J d'xF i 



2 \ i_-^qo „—ET 



(4.5) 



where I have used that at large (euclidean) time T the partition function picks up the 
ground state or vacuum energy E. For the sake of brevity I do not write the gauge 
fixing and Faddeev-Popov ghost terms. If the state is homogeneous, the energy can 
be written as E = O^V^ where 9^ v is the stress-energy tensor and is the 3- 
volume of the system. Hence, at large 4-volumes V = V^'T the partition function 
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is Z = exp(— 644V). This #44 includes zero-point oscillations and diverges badly. 
A more reasonable quantity is the partition function, normalized to the partition 
function understood as a perturbative expansion about the zero-field vacuum^, 

-9™-)V]. (4.6) 



exp 



Zp.T. 

We expect that the non-perturbative vacuum energy density 6*44— 9^ T ' is a negative 
quantity since we have allowed for tunneling: as usual in quantum mechanics, it lowers 
the ground state energy. If the vacuum is isotropical, one has #44 = 6*^/4. Using the 
trace anomaly, 



/V) fe) 2 ^ _ h I? 
4# 4 32tt 2 32tt 

where (3(g 2 ) is the Gell-Mann-Low function 



9 = Muy ; v ^ s b-^- (A 7) 



one gets 



2 exp ( b -V(Fl/32n 2 ) NP ) (4.9) 



Zp.T. V 4 

where {F 2 u )^p is the gluon field vacuum expectation value which is due to non- 
perturbative fluctuations, i.e. the gluon condensate. The aim of any QCD-vacuum 
builder is to minimize the vacuum energy or, equivalently, to maximize the gluon 
condensate. 

It is important that it is a renormalization-invariant quantity f], meaning that its 
dependence on the ultraviolet cutoff M and the bare charge g 2 {M) given at this cutoff 
is such that it is actually cutoff- independent: 



2 The latter can be distinguished from the former by imposing a condition that it does not contain 
integration over singular Yang-Mills potentials; recall that the instanton potentials are singular at 
the origins. 

3 To be more precise, the renorm-invariant quantity is (Ouu), see eq. (4.7); however if the coupling 
constant at the ultra-violet cutoff scale g 2 (M) is small enough, it is sufficient to use the beta function 
from one loop. 
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(F 2 J32tt 2 ) np = c Mexp(-£ 



r9 2 (M) dg 2 

w 



1 '1 



c'M exp 



32vr 2 



(4.10) 



bg 2 (M) 

By definition of the QCD partition function ( |4.5|) , the l.h.s. of eq. ( |4.10| ) is equal 
to -^d\nZ/d[87r 2 /g 2 (M)}. Applying the same differentiation operator to eq. (|4 . 1 0|) 
one gets a low-energy theorem p2]| : 



d 2 In (Z/Zp.t.] 



(d 



8tt 2 



g 2 (M) 



F 2 

d A x^- 



327T 2 



F 2 



32tt 2 



F 2 
d 4 x- ^ 



32tt 2 



4 
b 



d A x 



F 2 
32^ 



(4.11) 



If the bare coupling g 2 (M) is not chosen small enough, there are obvious corrections 
to this formula, following from the higher-order terms in the beta function. 

This low-energy theorem has an instructive consequence for instantons, predicting 
the dispersion of the number of pseudoparticles in a given 4-dimensional volume 0. 
Assuming the instanton ensemble is sufficiently dilute (corrections to this assumption 
will be discussed below) one can rewrite the low-energy theorem (|4.11[ 



as 



(N 2 ) - {N)< 



■(N) 
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(4.12) 



b" ' l\N c 

where = N + + N_ is the total number of Fs and Fs . 

Thus it follows from the renormalization properties of the Yang-Mills theory that 
the dispersion of the number of pseudoparticles is less than for a free gas for which one 
would get a Poisson distribution with (A 2 ) — (A) 2 = (A). In particular, at A c — ► oo 
the dispersion becomes zero, as it should be. 

One concludes that some kind of interaction of instantons with each other is crucial 
to support the needed renormalization properties of the underlying theory: any cutoff 
of the integrals over instanton sizes "by itself leads to the Poisson distrubution and 
hence to the violation of the low energy theorem ( [4.12|) . 

Differentiating In Z many times in respect to the bare coupling g 2 (M), one can 
easily generalize eq. ( |4.12j ) to any moments of the distribution. In short, the distru- 
bution in the number of pseudoparticles should be (for large (A)) 



P(N) ~ exp 



4 I (A) 



1 



(4.13) 
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5 One-instanton weight 



The words "instanton vacuum" mean that one assumes that the QCD partition func- 
tion is mainly saturated by an ensemble of interacting J's and J's together with quan- 
tum fluctuations about them. Instantons are necessarily present in the QCD vacuum 
if only because they lower the vacuum energy in respect to the purely perturbative 
(divergent) one. The question is whether they give the dominant contribution to the 
gluon condensate, and to other basic quantities. To answer this question one has to 
compute the partition function ( |4.5|) assuming that it is mainly saturated by instan- 
tons, and to compare the obtained gluon condensate with the phenomenological one. 
This work has been done a decade ago in ref. J7|; today direct lattice measurements 
confirm that the answer to the question is positive: the observed density of J's and 



J's is in agreement with the estimate ([4.4|) 



The starting point of this calculation is the contribution of one isolated instanton 
to the partition function ( |4.5|) or the one-instanton weight. It has been computed by 
't Hooft P|, and generalized to arbitrary groups by Bernard in a very clearly written 
paper |33| . 



The general field can be decomposed as a sum of a classical field of an instanton 
A^x, £) where £ is a set of 4iV c collective coordinates characterizing a given instanton 
(see eq. ( |3.13| )), and of a presumably small quantum field a M (x): 

A^x) = A I ll {x,£) + a li {x). (5.1) 

There is a subtlety in this decomposition due to the gauge freedom: an interested 
reader is addressed to ref. where this subtlety is treated in detail. 
The action is 

1 r 8ir 2 1 r 

Acti ° n = V J dx F i = + ~2 J d ^ x D » F ^ a » 

Here the term linear in drops out because the instanton field satisfies the 
equation of motion. The quadratic form W^ u has 4N C zero modes related to the 
fact that the action does not depend on 4N C collective coordinates. This brings in a 
divergence in the functional integral over the quantum field which, however, can 
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and should be qualified as integrals over the collective coordinates: centre, size and 
orientations. Formally the functional integral over a M gives 

, 1 (5.3) 
/det W^AT) 

which must be i) normalized (to the determinant of the free quadratic form, i.e. with 
no background field) , ii) regularized (for example by using the Pauli-Villars method) , 
and in) accounted for the zero modes. Actually one has to compute a "quadrupole" 
combination, 



det'W det(W + M 2 )l 5 

det W det{W + M 2 )\ ' ( ' ' 

where Wo is the quadratic form with no background field and M 2 is the Pauli-Villars 
mass playing the role of the ultra-violet cutoff; the prime reminds that the zero modes 
should be removed and treated separately. The resulting one-instanton contribution 
to the partition function (normalized to the free one) is 0, |33f : 



-inst 



d {p) 



Zp.T. 

p 5 



d A z^ J dp J dU d (p), 



57T 



-I 2N C 



9 2 (M) 



(M P y 



-N c 



8tt 2 



exp 



(5.5) 



g\M))- 

The product of the last two factors is actually a combination of the cutoff M and 
the bare coupling constant g 2 {M) given at this cutoff, which is cutoff-independent; 
it can be replaced by (Aq C dp) UNc ^ 3 , see eq. fl4.10|) . This is the way Aq C d enters 
into the game; henceforth all dimensional quantities will be expressed through Aqcd, 
which is, of course, a welcome message. The numerical coefficient C(N C ) depends 
explicitly on the number of colours; it also implicitly depends on the regularization 
scheme used. In the Pauli-Villars scheme exploited above |33| 



C{N C 



4.60exp(-1.68iV c 



vr 2 (iV c -i)!(iV c -2)f 
If the scheme is changed, one has to change the coefficient C(N C ) —>■ C'(N C ) 
(A/A') 6 . One has: A P . V . = lWA^s = 31.32A Zot = ... 

We have thus obtained the instanton weight in the one-loop approximation: 



(5.6) 
C(N C )- 
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d (p) = ( ^P(M) 2Nc expi-P(p)} 



(5.7) 



where (3(p) = 8n 2 /g 2 (p) is the one-loop inverse charge (not to be confused with the 
Gell-Mann-Low function!) 



000 =6 In 



1 



11 



(51 



M) 3 

Note that the /3 factor in the pre-exponent starts to "run" only at the 2- loop level, 
hence its argument is taken at the ultra-violet cutoff M. 

In the 2-loop approximation the instanton weight is given by |34| 



d (p) = ^(3(p) 2N ^x P 



-3^p) + (2N - h -\-^ 



0(1/ P{p)) 



, (5.9) 



where /3 n (p) is the inverse charge to the two-loop accuracy: 



^( P )=f3(p) + -\n 



3 c 



(5.10) 



Notice that both one- and two- loop eqs. (577, |5.9| ) formulae show that the integral 
over the instanton sizes p in eq. ( |5.5|) diverges as a high power of p at large p: this is 
of course the consequence of asymptotic freedom. It means that individual instantons 
tend to swell. This circumstance plagued the instanton calculus for many years. If one 
attemts to cut the p integrals "by hand", one violates the renormalization properties 
of the YM theory, as explained in the previous section. Actually the size integrals 
appear to be cut from above due to instanton interactions. 



6 Instanton ensemble 

To get a volume effect from instantons one needs to consider an II ensemble, with 
their total number N proportional to the 4-dimensional volume V. Immediately 
a mathematical difficulty arises: any superposition of J's and J's is not, strictly 
speaking, a solution of the equation of motion, therefore, one cannot directly use the 
semiclassical approach of the previous section. There are two ways two overcome this 
difficulty. One is to use a variational principle [0], the other is to use the effective 
YM lagrangian in the instanton field . 
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The idea of the variational principle is to use a modified YM action for which a 
chosen II ansatz is a saddle point. Exploiting the convexity of the exponent one can 
prove that the true vacuum energy is less than that obtained from the modified action. 
One can therefore use variational parameters (or even functions) to get a best upper 
bound for the vacuum energy. We call it the Feynman variational principle since the 
method was suggested by Feynman in his famous study of the polaron problem. The 
gauge theory is more difficult, though: one has not to loose either gauge invariance 
or the renormalization properties of the YM theory. These difficulties were overcome 
in ref. [0. A decade later I still do not think one can do an analytical evaluation of 
the II ensemble much better than in that paper: after all we are dealing with the 
"strong interactions", meaning that all dimenionless quantities are generally speaking 
of the order of unity - there are no small parameters in the theory. Therefore, one 
has to use certain numerical methods, and the variational principle is among the best. 
Todays direct lattice investigation of the II ensemble seem to indicate that Petrov 
and I have obtained rather accurate numbers in this terrible problem. 

The normalized (to perturbative) and regularized YM partition function takes the 
form of a partition function for a grand canonical ensemble of interacting psuedopar- 
ticles of two kind, Ps and Ps : 



collective coordinates of (anti)instantons: their coordinates z, sizes p and orientations 
given by SU(N C ) unitary matrices in the adjoint representation 0\ dO means the Haar 
measure normalized to unity. The instanton interaction potential Pint (to be discussed 
below) depends on the separation between pseudoparticles, z m — z n , their sizes p m) „ 
and their relative orientations O m O^. 

In the variational approach the interaction between instantons arise from i) the 
defect of the classical action, ii) the non-factorization of quantum determinants and 
Hi) the non-factorization of jacobians when one passes to integartion over the col- 
lective coordinates. All three factors are ansatz-dependent, but there is a tendency 
towards a cancellation of the ansatz-dependent pieces. Qualitatively, in any ansatz 
the interactions between Ps and Ps resemble those of molecules: at large separations 
there is an attraction, at smaller separations there is a repulsion. It is very important 
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that the interactions depend on the relative orientations of instantons: if one averages 
over orientations (which is the natural thing to do if the II medium is in a disordered 
phase; if not, one would expect a spontaneous breaking of both Lorentz and colour 
symmetries |7J), the interactions seem to be repulsive at any separations. 

In general, the mere notion of the instanton interactions is notorious for being 
ill-defined since instanton + antiinstanton is not a solution of the equation of motion. 
Such a configuration belongs to a sector with topological charge zero, thus it seems to 
be impossible to distinguish it from what is encountered in perturbation theory. The 
variational approach uses brute force in dealing with the problem, and the results 
appear to be somewhat dependent on the ansatz used. Thanks to the inequality for 
the vacuum energy mentioned above, we still get quite a useful information. However, 
recently a mathematically unequivocal definition of the instanton interaction poten- 
tial has been suggested, based on analyticity and unitarity |3(| ^5[. This definition 
automatically cuts off the contribution of the perturbation theory. The first three 
leading terms for the interaction potential at large separations has been computed 



[ |35|j , at smaller separations one observes a strong repulsion [37], though the exact 
form is still unknown. 

Summing up the discussion, I would say that today there exists no evidence that a 
variational calculation with the simplest sum ansatz used in ref. |7J is qualitatively of 
even quantitatively incorrect, therefore I will cite the numerics from those cacluations 
in what follows. 

The main finding is that the II ensemble stabilizes at a certain density related to 
the Aqcd parameter (there is no other dimensional quantity in the theory!): 



~ -(<%) -'-> (0.75A fis ) 4 (6.2) 



—{QD * H 

yV*TI y 

which would require A^g — 265 MeV to get the phenomenological value of the 
condensate. It should be mentioned however that using more sophisticated variational 
Ansatze one can obtain a larger coefficient in eq. ( |6.2| ) and hence would need smaller 
values of A. 

The average sizes p appear to be much less than the average separation R. Nu- 
merically we have found for the SU(3) colour: 

14 <-> 
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which coincides with what was suggested previously by Shuryak || from considering 
the phenomenological applications of the instanton vacuum. This value should be 
compared with that found from direct lattice measurements |T7| |: p/R ~ .37 — .4, 
depending on where one stops the cooling procedure. The packing fraction, i.e. the 
fraction of the 4-dimensional volume occupied by instantons apears thus to be rather 
small, 7r 2 p 4 /^ 4 ~ 1/8. This small number can be traced back to the "accidentally" 
large numbers appearing in the 4-dimensional YM theory: the lliV c /3 of the Gell- 
Mann-Low beta function and the number of zero modes being 4N C . We have checked 
that the same variational principle applied to the 2-dimensional sigma models also 
possessing instantons, does not yield a small packing fraction. The 4-dimensional 
YM theory seems to be simpler from this angle. Meanwhile, it is exactly this small 
packing fraction of the instanton vacuum which gives an a posteriori justification for 
the use of the semi-classical methods. As I shall show in the next sections, it also 
enables one to identify adequate degrees of freedom to describe the low-energy QCD. 



7 Chiral symmetry breaking in QCD 

The QCD lagrangian with Nf massless flavours is known to posses a large global 
symmetry, namely a symmetry under U(Nf) x U(Nf) independent rotations of left- 
and right-handed quark fields. This symmetry is called chiral. Instead of rotating 
separately the 2-component Weyl spinors one can make independent vector and axial 
U(Nf) rotations of the full 4-component Dirac spinors - the QCD lagrangian is in- 
variant under these transformations too. Meanwhile, the axial transformation mixes 
states with different P-parities. Therefore, if that symmetry remains exact one would 
observe parity degeneracy of all states with otherwise the same quantum numbers. 
In reality the splittings between states with the same quantum numbers but opposite 
parities is huge. For example, the splitting between the vector p and the axial a\ 
meson is (1200 — 770) ~ 400 MeV; the splitting between the nucleon and its parity 
partner is even larger: (1535 — 940) ~ 600 MeV. (Another possibility is that the 
nucleon is just massless, which looks even less pleasant.) 

The splittings are too large to be explained by the small bare or current quark 
masses which break the chiral symmetry from the beginning. Indeed, the current 
masses of light quarks are: m u ~ 4 MeV, ~ 7 MeV, m s ~ 150 MeV. The 
only conclusion one can draw from these numbers is that the chiral symmetry of the 
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QCD lagrangian is broken down spontaneously, and very strongly. Consequently, one 
should have light (pseudo) Goldstone pseudoscalar hadrons - their role is played by 
pions which indeed are by far the lightest hadrons. 

The order parameter associated with chiral symmetry breaking is the so-called 
chiral or quark condensate: 



(V>V) - -(250 MeVf. (7.1) 

It should be noted that this quantity is well defined only for massless quarks, otherwise 
it is somewhat ambigious. By definition, this is the quark Green function taken at one 
point; in momentum space it is a closed quark loop. If the quark propagator has only 
the "slash" term, the trace over the spinor indices implied in this loop would give an 
identical zero. Therefore, chiral symmetry breaking implies that a massless (or nearly 
massless) quark develops a non-zero dynamical mass (i.e. a "non-slash" term in the 
propagator). There are no reasons for this quantity to be a constant independent 
of the momentum; moreover, we understand that it should anyhow vanish at large 
momentum. The value of the dynamical mass at zero momentum can be estimated 
as one half of the p meson mass or one third of the nucleon mass, that is about 
M(0) ~ 350 — 400 MeV; this scale is also related to chiral symmetry breaking and 
should be emerge together with the condensate (|T7T| ) . 

One could imagine a world without confinement but with chiral symmetry break- 
ing: it would not be drastically different from what we meet in reality. There would 
be a tightly bound light Goldstone pion, and relatively loosely bound p meson and 
nucleon with correct masses, which, however, would be possible to ionize from time to 
time. Probably the spectrum of the highly excited hadrons would be wrong, though 
even that is not so clear |BS|. We see, thus, that the spontaneous chiral symmetry 
breaking is the main dynamical happening in QCD, which determines the face of the 
strong interactions world. In the next sections I explain why and how chiral symme- 
try is broken in the instanton vacuum, and why it is a most realistic picture. The 
forthcoming sections are based on our work with Petrov [H, HI. 
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8 Chiral symmetry breaking by definition 



I start by writing down the QCD partition function. Functional integrals are well 
defined in euclidean space which is obtained by the following formal substitutions of 
Minkowski space quantitites: 



i%Mo — xea, %Mi — %Ei, A M0 — iA EA , A Mi — A Ei , 

i$M = We, Imo = 7£4, 7m; = ijEi, 1m 5 = lE5- 

Neglecting for brevity the gauge fixing and Faddeev-Popov ghost terms 
partition function can be written as 

I" 1 Nf 
Z = J DA^Dtf exp -— J F% + £ J ^}(*V + im f )^ f 



(8.1) 
the QCD 



J DAf, exp 



1 



4g2 I ^ 



Nt 



Y[det(iV + im f ). 
J / 



The chiral condensate of a given flavour / is, by definition 

1 d 



The Dirac operator has the form 



V drrif 



(8.2) 



.3) 



iV = 7 ^ + 4 + « ( 8 - 4 ) 
where A 1 ^ denotes the classical field of the II ensemble and is a presumably small 
field of quantum fluctuations about that ensemble, which I shall neglect as it has 
little impact on chiral symmetry breaking. Integrating over DA^ in eq. ( |3.2| ) means 
averaging over the II ensemble, therefore one can write 



Z = dei{iV + im) (8.5) 

where I temporarily restrict the discussion to the case of only one flavour for simplicity. 
Because of the im term the Dirac operator in flS.5|) is formally not hermitean; however 
the determinant is real due to the following observation. Suppose we have found the 
eigenvalues and eigenf unctions of the Dirac operator, 
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(8.6) 



then for any A n 7^ there is an eigenfunction <3v = l^n whose eigenvalue is 
An' = — A n . This is because 75 anticommutes with iV. Owing to this the fermion 
determinant can be written as 



det(zV + im) = ]~J(A n + im) = y JJ(A^ + m 2 ) = exp 



exp 



£ j —00 



d\ u(X) ln(A 2 + m 2 



•7) 



where I have introduced the spectral density z/(A) of the Dirac operator zV. Note that 
the last expression is real and even in m, which is a manifestation of the QCD chiral 



mvanance. 



Differentiating eq. ( |8.7| ) in m and putting it to zero one gets according to the 
general eq. ( |3.3|) a formula for the chiral condensate: 



m- 1 



V dm 



- / d\ z/(A) ln(A 2 + m 2 ) = -jt d\v(\) 



in 



A 2 + m 2 ) 



m->0 

(8-8) 

where z/(A) means averaging over the instanton ensemble together with the weight 
given by the fermion determinant itself. The latter, however, may be cancelled in the 
so-called quenched approximation where the back influence of quarks on the dynamics 
is neglected. Theoretically, this is justified at large N c . 

Naively, one would think that the r.h.s. of eq. ( |B.8|) iz zero at m — > 0. That would 
be correct for a finite-volume system with a discrete spectral density. However, if the 
volume goes to infinity faster than m goes to zero (which is what one should assume 
in the thermodynamic limit) one should use instead 

"' '^ >0 sign(m)7T(5(A), (8.9) 



A 2 + m 
so that one gets p9| : 



(ifnj)) = __ s ign(m)7rz/(0). (8.10) 

The chiral condensate is thus proportional to the averaged spectral density of 
the Dirac operator at zero eigenvalues. The appearance of the sign function is not 
accidental: it means that at small m QCD partition functions depends on m non- 
analytically: 
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\nZ = V(cq + 7cu(0)\m\ + C2m 2 ln(|m|) + 



The fact that the partition function is even in m is the reflection of the original 
invariance of the QCD under 75 rotations; the fact that it is non-analytic in the 
symmetry-breaking parameter m is typical in the situation where symmetry is broken 
spontaneously 

A generalization of the above formulae to the case of several flavours is simple 
|tj. Consider Nf quark flavours with a most general mass matrix 



. 1+75 , . 1 
%m L — V im,R— 



75 



(8.12) 



The fermion determinant can be written as 
det(zV + im) = exi 



exp 



■ / dXu(X) In det^ (A 2 + m^m^ 
K^) J -7 T^Nf exp (-£(A 2 + m L m R 



(8.13) 



where A 2 is a Nf x Nf matrix proportional to the unit matrix. Let us expand i/(A) at 
small A, and leave only the z/(0) term. Integrating ( P-13| ) in A we get the correspondent 
non-analytic term in the partition function: 



/qr poo Hf 

\ nZ ~ -!L-v(0)J -^TT Nf [exp(-tm L m R )} = nu(0)Ti Nf (m L m R ) 1 / 2 . (8.14) 

This expression is non-analytic in the mass matrix; differentiating it in respect to the 
masses one finds the phases of the condensate. 

In the next sections I shall present three different derivations of the fact that 
instantons indeed break chiral symmetry. 

9 Chiral symmetry breaking by instantons: quali- 
tative derivation 

The key observation is that the Dirac operator in the background field of one (anti) 
instanton has an exact zero mode with A = @. It is a consequence of the general 
Atiah-Singer index theorem; in our case it is guaranteed by the unit Pontryagin index 
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or the topological charge of the instanton field. These zero modes are 2-component 
Weyl spinors: right-handed for instantons and left-handed for antiinstantons. Explic- 
itly, the zero modes are (a = l...N c is the colour and i,j, k = 1,2 are the spinor 
indices): 

- z 2 )\ a i = 4>(x - z 2 , p 2 )(x - z 2 )r j U^ k e jk , 

P) = 7r(2a; 2 ) 1 /2( x 2 + p 2 ) 3/2- ( 9 - X ) 
Here z^, P\,U\ are the center, size and orientation of an instanton and z 2fl , p 2 , U 2 are 
those of an antiinstanton, respectively, e jfc is the 2x2 antisymmetric matrix. 

For infinitely separated I and I one has thus two degenerate states with exactly 
zero eigenvalues. As usual in quantum mechanics, this degeneracy is lifted through 
the diagonalization of the hamiltonian, in this case the hamiltonian is the full Dirac 
operator. The two "wave functions" which diagonalize the "hamiltonian" are the sum 
and the difference of the would-be zero modes, one of which is a 2-component left- 
handed spinor, and the other is a 2-component right-handed spinor. The resulting 
wave functions are 4-component Dirac spinors; one can be obtained fron another by 
multiplying by the 75 matrix. As the result the two would-be zero eigenstates are split 
symmetrically into two 4-component Dirac states with non-zero eigenvalues equal to 
the overlap integral between the original states: 

A = ±|T 12 |, T 12 = [ d 4 x<S>{(-i jd)® 2 Rl ^T -^lTr(UlU 2 Rt 2 ). (9.2) 
J K 12 

We see that the splitting between the would-be zero modes fall off as the third 
power of the distance between I and J; it also depends on their relative orientation. 
The fact that two levels have eigenvalues ±A is in perfect agreement with the 75 
invariance mentioned in the previous section. 

When one adds more J's and J's each of them brings in a would-be zero mode. 
After the diagonalization they get split symmetrically in respect to the A = axis. 
Eventually, for an iT ensemble one gets a continuous band spectrum with a spectral 
density u(X) which is even in A and finite at A = 0. 

One can make a quick estimate of v(0): Let the total number of J's and J's in the 
4-dimensional volume V be N. The spread A of the band spectrum of the would-be 
zero modes is given by their average overlap ( |9.2|) : 
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&~^j{tfRlV)T{R)T*{R)~-L (9.3) 

where p is the average size and R = (N/V)~ l l A is the average separation of instantons. 
Note that the spread of the would-be zero modes is parametrically much less than l/p 
which is the typical scale for the non-zero modes. Therefore, neglecting the influence 
of the non-zero modes is justified if the packing fraction of instantons is small enough. 
From eq. ( 8.10|) one gets an estimate for the chiral condensate induced by instantons: 



<#) = -£,(0) - -~ ~ (9.4) 

Note that the chiral condensate appears to be proportional to the square root of 
the instanton density N/V: again it is as it should be for the order parameter of 
spontaneous symmetry breaking. 

It is amusing that the physics of the spontaneous breaking of chiral symmetry 
resembles the so-called Mott-Anderson conductivity in disordered systems. Imagine 
random impurities (atoms) spread over a sample with final density, such that each 
atom has a localized bound state for an electron. Due to the overlap of those localized 
electron states belonging to individual atoms, the levels are split into a band, and the 
electrons become delocalized. That means conductivity of the sample. In our case the 
localized zero quark modes of individual instantons randomly spread over the volume 
get delocalized due to their overlap, which means chiral symmetry breaking. 

There is a difference with the Mott-Anderson conductivity, though. In the case of 
atoms, the wave functions have an exponential falloff, so that the overlap integrals are 
exponentially small at large separations. It means that, if the density of impurities is 
small enough, there might be a phase transition to an insulator state. In our case the 
eigenvalue is exactly zero, and the zero modes decay as a power - see eq. Q9.1|) . The 
overlap integrals are also decreasing as a power of the separation (see eq. (|9.2p ), and 
chiral symmetry breaking occurs at any instanton densities - at least in the quenched 
approximation, when one neglects the back influence of quarks on the dynamics of 
instantons. However, at non-zero temperatures the density of instantons decrease 



and the fermion zero modes have an exponential falloff therefore one can expect 
a Mott-Anderson phase transition to an insulator state, meaning the restoration of 
chiral symmetry. It should be noted that chiral symmetry may be restored as due to 
the back influence of quarks on the instanton ensemble - it is an 0(Nf/N c ) effect then. 
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It has been studied recently in refs. [[43]] where a formation of instanton molecules has 
been suggested as a mechanism of chiral symmetry restoration at high temperatures. 
In fact, the restoration mechanism may be different for different numbers of colours 
and flavours. 

I should mention that the idea that instantons can break chiral symmetry has 
been discussed previously (see refs. [44, [45], [5], 0) however the present mechanism 
and a consistent formalism has been suggested and developed in papers || || . 

Recently there have been much interesting work done generalizing this mecha- 
nism of chiral symmetry breaking to other configurations [[4(| and studying general 
properties of the spectral density of the Dirac operator for various ensembles ETl [48 



10 Derivation II: quark propagator 

Having explained the physical mechanism of chiral symmetry breaking as due to the 
derealization of the would-be zero fermion modes in the field of individual instantons, 
I shall indicate how to compute observables in the instanton vacuum. The main 
quantity is the quark propagator in the instanton vacuum, averaged over the instanton 
ensemble. This quantity has been calculated in refs. |8|, [HJ. In particular, Pobylitsa 
has derived a closed equation for the averaged quark propagator, which can be 
solved as a series expansion in the formal parameter Np A /VN c which numerically is 
something like 1/250. 



The result of refs. 0, [K| is that in the leading order in the above parameter the 
quark propagator has the form of a massive propagator with a momentum- dependent 
dynamical mass: 



where F(z) is a combination of the modified Bessel functions and is related to the 
Fourier transform of the zero mode ( |SrJ. 1[ ) : it is equal to one at z = and decreases 
rapidly with the momentum, measured in units of the inverse average size of instan- 
tons (see next section); c is a constant of the order of unity which depends slightly on 
the approximation used in deriving the propagator. Note that the dynamical quark 
mass is non-analytical in the instanton density (similar to the chiral condensate it is 
an order parameter for spontaneous symmetry breaking). 
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Fixing the average density by the empirical gluon condensate (see section 4) so 
that R ~ 1 fm and fixing the ratio p/R = 1/3 from our variational estimate, we get 
the value of the dynamical mass at zero momentum, 

M(0) ~ 350 MeV (10.2) 

while the quark condensate is 



m = i I gLvSM = -4N C J ^j^fej - -(255 MeV f. (10.3) 



Both numbers, ( 10.2j ) and ( |10.3| ), appear to be close to their phenomenological 
values. 

Using the above small parameter one can also compute more complicated quanti- 
ties like 2- or 3-point mesonic correlation functions of the type 

(J A (x)J B (y)), (J A (x)J B (y)J c (z)), J A = $T A 1> (10.4) 

where T A is a unit matrix in colour but an arbitrary matrix in flavour and spin. 
Instantons influence the correlation function in two ways: i) the quark and antiquark 
propagators get dressed and obtain the dynamical mass, as in eq. ( |10.1| ), ii) quark 
and antiquark may scatter simultaneously on the same pseudoparticle; that leads to 
certain effective quark interactions. These interactions are strongly dependent on the 
quark- antiquark quantum numbers: they are strong and attractive in the scalar and 
especially in the pseudoscalar and the axial channels, and rather weak in the vector 
and tensor channels. I shall derive these interactions in the next section, but already 
now we can discuss the pseudoscalar and the axial isovector channels. These are the 
channels where the pion shows up as an intermediate state. 

Since we have already obtained chiral symmetry breaking by studying a single 
quark propagator in the instanton vacuum, we are doomed to have a massless Gold- 
stone pion in the appropriate correlation functions. However, it is instructive to follow 
how does the Goldstone theorem manifest itself in the instanton vacuum. It appears 
that technologically it follows from a kind of detailed balance in the pseudoscalar 
channel (such kind of equations are encountered in perturbative QCD where there is 
a delicate cancellation between real and virtual gluon emission). However, since we 
have a concrete dynamical realization of chiral symmetry breaking we can not only 
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check the general Ward identities of the PCAC (which work of course) but we are in 
a position to find quantities whose values do not follow from general relations. One of 
the most important quantities is the F v constant: it can be calculated as the residue 
of the pion pole. We get: 



const / p 



P 



7i\ 2 



RJ \ 



R 

In— ~100MeV vs. 93 MeV (exper.). (10.5) 
P 



This is a very instructive formula. The point is, F n is surprisingly small in the strong 
interactions scale which, in the instanton vacuum, is given by the average size of 
pseudoparticles, 1/p — 600 MeV. The above formula says that F^ is down by the 
packing fraction factor (p/R) 2 ~ 1/9. It can be said that F n measures the diluteness 
of the instanton vacuum! However it would be wrong to say that instantons are in a 
dilute gas phase - the interactions are crucial to stabilize the medium and to support 
the known renormalization properties of the theory, therefore they are rather in a 
liquid phase, however dilute it may turn to be. 

By calculating three-point correlation functions in the instanton vacuum we are 
able to determine, e.g. the charge radius of the Goldstone excitation: 

r 2 ^ v^£ ^ (340 MeV)' 1 vs. (310 MeV)- 1 (exper.). (10.6) 
2-kFk; 

Let me note that all quantities exhibit the natural behaviour in the number of 
colours N c : 

(*%>)" y = 0{N e ), (M = 0(N e ), F 2 = 0(N c ) } 

p = 0(l), M(0) = O(l), V^ = 0(1), etc. (10.7) 
A systematic numerical study of various correlation functions in the instanton 



vacuum has been performed by Shuryak, Verbaarschot and Schaefer |fL5|| , see also 
Shuryak's lectures at this School. In all cases considered the results agree well or very 
well with experiments and phenomenology. As I already mentioned in the introduc- 
tion, similar conclusions have been recently obtained from direct lattice measurements 



[fL7|j . I think that one can conclude that instantons are explaining the basic properties 



of the QCD ground state and that of light hadrons. 
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11 Derivation III: Nambu— Jona-Lasinio model 



The idea of the first two derivations of chiral symmetry breaking by instantons, pre- 
sented above, is: " Calculate quark observables in a given background gluon field, then 
average over the ensemble of fields", in our case the ensemble of J's and J's . The 
idea of the third derivation is the opposite: "First average over the II ensemble and 
obtain an effective theory written in terms of interacting quarks only Then compute 
observables from this effective theory". This approach is in a sense more economical; 
it has been developed in refs. || £4(J . 

Quark interaction arises when two or more (anti) quarks happen to scatter over 
the same pseudoparticle; averaging over its positions and orientations results in a 
four- (or more) fermion interaction term whose range is that of the average size of 
instantons. The most essential way how instantons influence quarks is, of course, 
via the zero modes. Since each massless quark flavour has its own zero mode, it 
means that the effective quark interactions will be actually 2Nf fermion ones. They 
are usually referred to as 't Hooft interactions as he was the first to point out the 
quantum numbers of these effective instanton- induced interactions f2|. In case of 
two flavours they are four-fermion interactions, and the resulting low-energy theory 



resembles the old Vaks-Larkin-Nambu- Jona-Lasinio model ||9| which is known to 
lead to chiral symmetry breaking. In this section I derive this model from instantons, 
following refs. |9|, |4(J . Recently it has been revisited in ref . |50 . 



The starting point is the quark Green function in the field of one instanton. It 
can be written as a sum over all eigenf unctions of the correspondent Dirac operator 
& n (x) - see eq. flSlf): 

V K + im 

= -*°w«K») +gf( } (1L1) 

im 

where $ is the zero mode ( |9.1[ ) (henceforth I omit the subscript 0), and S' is 
the sum over non-zero modes, which is finite in the chiral limit m — > 0. For the 
simplicity of the derivation we replace it by the free Green function So(x,y), though 
the exact propagator in the field of one instanton is known. Thus instead of the exact 
propagator we write 
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S'(x, y) = - $(x)$t(y) + S (x, y). (11.2) 
im 

This approximate Green function is correctly taking into account the zero mode, 
that is the low-momentum part, and at large momentum it reduces to the free Green 
function, as it should. Therefore, it is an interpolation of the exact propagator; at mo- 
menta p ~ 1/p the numerics will be not exact. However, phenomena related to chiral 
symmetry breaking correspond to lower momenta, and the use of the simplification 



11. 2|) is therefore theoretically justified. 

We have now to build the Green function (and other quantities) in the field of N + 
J's and AL J's and to average over their ensemble. To that end we use the following 
mathematical trick. Consider a fermion action 

exp {-A 1 ^,^}) = exp ^ d x ipH JQi^j (im - V 1 ^ ,tp]) , 

V 1 ^, V] = / d * x jd&ix)) J d A y (V%)z jdip(y)) (11.3) 

where $ 7 is the zero mode in the field of the Jth pseudoparticle. In case of Nf > 1 
one has to take here a product of the {irrif — V 1 ) factors for all flavours. 

The action ( |11.3|) has the following properties: i) the correspondent partition 



function normalized to the free one is equal to im, as it should be in case of one 
instanton; ii) the Green function, computed with this action coincides with that of 



eq. (|11.2|) . It means that this action correctly describes quarks in the field of a given 
instanton at low and at large momenta, and interpolates in between. 
In the field of iV+ J's and iV_ J's the fermion action is 



/ r \ N++N - 

exp(-A[^]) = rjexp I tfxrfiMf) 1] [imf-V 1 ^,^]). (11.4) 

To get the QCD partition function one has to integrate over the quark fields and 
average over the ensemble of instantons: 

Zqcd = J D^D^(exp(-A[iij\^])) (11.5) 

where (...) denotes averaging over the ensemble. It was shown in ref. [0] (see also 
[ |50|j ) that the J J ensemble can be described by an effective one-particle distribution 
in the instanton sizes, which can be found from a variational principle: 
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d(p) = const p b 5 exp 



p 2 b-4: 



3 3 ; 



11.6) 



r 2 ; ' 

At large N c it is 5-peaked around the average p. Therefore, one can replace the 
averaging over ensemble by independent averaging over positions and orientations of 
invidual pseudoparticles. Thus eq. ( |11.5| ) becomes 



Zqcd = J DipDip^exv (J d A xip\i /% 



11.7) 



■ n [im f - V[$, ifrf] J (im f - V[$, V/] 

where the bar means averaging over individual pseudoparticles. For simplicity I shall 
replace all sizes by their average p. I shall also consider only the chiral limit rrif — > 0. 
It is natural to introduce non-local 2Nf fermion vertices: 



Y + 



J f 



11.8) 



where V 1 ^ depends on the (anti) instanton centers z M and orientations U through 
the zero modes see eq. ( |11.3| ) and ( |9.1| ). The partition function ( |11.7| ) can be 

written as 



Zqcd = J DipfDipjexj) f^J d 4 x^ f i /)ip f 



d\^ 



2tt 



J dT± exp 



iV+ln-^ + 



(11.9) 



Indeed, integarting over A± one gets 8(Y—F), and after integrating over T one recovers 
eq. ( |11.7| ). In the thermodynamic limit N±, V — > oo, N/V fixed, integration over 11 
and A± can be performed by the saddle point method. We integrate first over II: 



Z, 



QCD 



d\± 



exp 



N. + In 



i\+V 



-1 + + 



-) 



• J LhjjfDtj)) exp (J d 4 xip}i jd^f + i\+Y+ + i\_Y^j . (11.10) 

Because of the non-locality it is more convenient to write the vertices ( |11.8| ) in 
the momentum space. Let us decompose the 4-component Dirac spinors describing 
quark fields into left- and right-handed 2-component Weyl spinors which we denote 



as 
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^ukv ^HR)fai> ( 1L11 ) 

where / = l...Nf stand for flavour, a = l...N c stand for colour and i — 1,2 stand 
for spin indices. Let us introduce the formfactor functions F(k) which are related to 
the Fourier transforms of the zero modes (|9.1| ) and are attributed to each fermion 
entering the vertex (z = kp/2): 

F{k) = 2z[I (z)K 1 (z) - h(z)K (z) - h x {z)Ki{z) ^ J^, F(0) = 1. (11.12) 
The Nf fermion vertices can be written as: 

Y$, = j dU \[j $^ 2n P F ^ J ^^pF(ln)(2nr5(k 1 + ... + k Nf -h-...-l Nf ) 

■ u« : u^e^e 3nPn {-^ Lfnanin {KW L n Snjn (Q) ; (11.13) 

for the vertices induced by J's one has to replace left-handed Weyl spinors by 
rigth-handed ones. The integral J dU means averaging over instanton orientations in 
colour space. In particlular, one has 

JdU = l, JdUU*Ur = PK etc. (11.14) 

To get the 2Nf vertices in a closed form one has to perform explicitly integration 
over the instanton orientations. I present below the results U for Nf = 1,2 and for 
any Nf but iV c — > oo. 

Nf=± 



In this case the "vertex" Q11.8 ) is just a mass term: 



Yf = J -0^^(k)(2npF(k)r^^m- (H.15) 



One has to plug it into the eq. (|11.10|) , integrate over fermions, and find the saddle- 



point values of \±. If the CP symmetry is conserved, one has N + = N_ = N/2, and 
the saddle-point values satisfy A + = A_. Then the 75 term in Y ± get cancelled, and 



eq. ( 11.15Q gives a momentum-dependent mass 

M{k) = M(0)F 2 (k) } M(0) = X(2tt P ) 2 , (11.16) 
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where M(0) or A is found from the equation ||, |9] (called sometimes self-consistency 
or gap equation): 



AN C r d*k M 2 {k) 



= 1. (11.17) 

N/V J (2tt) 4 k 2 + M 2 {k) { } 

Let me mention that exactly the same gap equation ( |11.17] ) has been obtained M 
in another approach: by first finding the quark propagator in the instanton vacuum 
and then averaging over the instanton ensemble. It is seen from eq. ( |1 1 . 1 7| ) that 
the dynamically generated mass is of the order of M(0) ~ \J N/ (VN C ) p. Knowing 
the form of M(k) given by eq. ( |11.12j ) and using the "standard" values N/V = 
(1 fm)~ 4 , p = (1/3) fm we find numerically M(0) ~ 350 MeV. If one neglects M 2 



in the denominator of eq. ( |1 1 . 1 7| ) flOf one gets M(0) ~ 420 MeV. This deviation 
indicates the accuracy of the "zero mode approximation" used in this derivation: it 
is about 15%. 



N 



f 



In this case averaging over the instanton orientations gives a nontrivial 4-fermion 
interaction: 

2N 2 [d^k x d^k 2 d%d% 



Y, 



(+) 



N/V 



chh 



2(N 2 - 1 
1 



2N 



(2tt) 12 
- 1 



M(k x )M{k 2 )M{l x )M{h 



2N 



(^l /l (A: 1 )^Vf(Zi))(4/ 2 (^)^^f(Z2)) 



11.18) 



For the antiinstanton-induced vertex ^ one has to replace left handed components 



by right-handed. In eq. ( |11.18| ) I have included the factor A + and, moreover, fixed it 
from the saddle-point equation. With the normalization of eq. ( |11.18| ) the dynamical 
mass M{k) satisfies exactly the same gap equation ( |11.17| ) as in the case of Nf = 1. 
Note that the second (tensor) term is negligible at large N c . Using the identity 



2e hh 



'9192 



5 h 5 h 

91 92 



one can rewrite the leading (first) term of eq. (|11.18|) as 



11.19) 



11.20) 
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which resembles closely the Nambu-Jona-Lasinio model. It should be stressed though 
that in contrast to that at hoc model the interaction ( |11.18| ) i) violates explicitly the 
Ua{1) symmetry, ii) has a fixed interaction strength and Hi) contains an intrinsic 
ultraviolet cutoff due to the formfactor function M(k). This model is known to lead 
to chiral symmetry breaking, at least at large N c when the use of the mean field 
approximation to the model is theoretically justified. 

Any Nt 



At arbitrary Nf the leading term at N c — > oo can be written as a determinant of 
a Nf x Nf matrices composed of quark bilinears: 

J S )(X) = / ^^-^ v /M(A;)M(/)^(fc)^^(0- (H.21) 

Again, one can prove that at least at large N c , this interaction leads to the spontaneous 
chiral symmetry breaking, with the dynamical mass determined by the gap equation 
Q11.17 ), and the chiral condensate given by eq. (|10.3|) . The bosonization of these 



interactions has been performed in ref. ||; it paves the way to studying analytically 
various correlation functions in the instanton vacuum. 

A separate issue is the application of these ideas to hadrons made of heavy quarks 



12] and of light and heavy ones [51 



12 QCD at still lower energies 

Using the packing fraction of instantons p/R ~ 1/3 as a new algebraic parameter one 
observes that all degrees of freedom in QCD can be divided into two categories: i) 
those with masses > 1/p and ii) those with masses < If one restricts oneself to 
low-energy strong interactions such that momenta are < 1/p ~ 600 MeV, one can 
neglect the former and concentrate on the latter. There are just two kind of degrees 
of freedom whose mass is much less than the inverse average size of instantons: the 
(pseudo) Goldstone pseudoscalar mesons and the quarks themselves which obtain 
a dynamically-generated mass M ~ (l/p)(p 2 /R 2 ) <C 1/p. Thus in the domain of 
momenta k <C 1/p QCD reduces to a remarkably simple though nontrivial theory of 
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massive quarks interacting with massless or nearly massless pions. It is given by the 
partition function || H 



2 low mom. = J D ^ D ^ exp J d i x ^\ ^ ij9 + iMe™^ 15 ) 1p 



;i2.i) 



Notice that there is no kinetic energy term for the pions, and that the theory is not 
a renormalizable one. The last circumstance is due to the fact that it is an effective 
low-energy theory; the ultraviolet cutoff is actually 1/p. 

There is a close analogy with solid state physics here. The microscopic theory of 
solid states is QED: it manages to break spontaneously the translational symmetry, 
so that a Goldstone excitation emerges, called the phonon. Electrons obtain a "dy- 
namical mass" m* due to hopping from one atom in a lattice to another. The "low 
energy" limit of solid state physics is described by interactions of dressed electrons 
with Goldstone phonons. These interactions are more or less fixed by symmetry con- 
siderations apart from a few constants which can be deduced from experiments or 
calculated approximately from the underlying QED. Little is left of the complicated 
dynamics at the atom scale. 

What Petrov and I have attempted, is a similar path: one starts with the funda- 
mental QCD, finds that instantons stabilize at a relatively low density and that they 
break chiral symmetry; what is left at low momenta are just the dynamically massive 
quarks and massless pions. One needs two scales to describe strong interactions at 
low momenta: the ultra-violet cutoff, whose role is played by the inverse instanton 
size, and the dynamical quark mass proportional to the square root of the instanton 
density. If one does not believe our variational calculations of these quantities one 
can take them from experiment. 

If one integrates off the quark fields in eq. ( |12.1 ) one gets the effective chiral 
lagrangian: 

S eff [ir] = -N c Tt \n(ijd + iMW 5 ) , 

U = exp^vrV 4 ), W 5 = exp(z7r A r A 7 5 ), L pu = %U^d p U. (12.2) 
One can expand eq. (|12.2 ) in powers of the derivatives of the pion field and get 



P 2 r IV r 

SeffM = -f J d A x TrLj - -^- 2 J d A x [2Tr(^L M ) 2 + TiL /1 L u L fl L 
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]\[ c r 

+ J d x 6 <xfa6e TrLaLpLjLsLe + ... (12.3) 

The first term here is the old Weinberg chiral lagrangian with 

f° = if , f^L M2 < fc » ■ (12 4) 

the second term are the four- derivative Gasser-Leutwyler terms (with coefficients 
which turn out to agree with those following from the analysis of the data); the last 
term in eq. ( 12.3Q is the so-called Wess-Zumino term. Note that the F n constant 



diverges logarithmically at large momenta; the integral is cut by the momentum- 
dependent mass at k ~ 1/p, so that one gets the same expression as in a different 
approach described in section 10, see eq. ( |10.5| ). 



An ideal field of application of the low-momentum partition function ( 12.1| ) is the 



quark-soliton model of nucleons ]TT] - actually the model has been derived from this 
partition function. The size of the nucleon ~ (250 MeV) -1 is much larger than the 
size of instantons ~ (600 MeV)~ l ; hence the low-momentum theory Q12.1 ) seems to 



be justified. Indeed, the computed static characteristics of baryons like formfactors, 
magnetic moments, etc., are in a good accordance with the data (for a review see 



ref. IT^l). What is not yet computed, are the nucleon structure functions at a low 
normalization point, however given the previous experience, I can bet it would go 
through the data points. 



13 How instantons may help confinement 

Our analytical calculations sketched in these lectures, the extensive numerical studies 
of the instanton vacuum by Shuryak and collaborators and the recent direct lattice 
measurements - all point out that instantons play a crucial role in determining the 
world of light hadrons, including the nucleon. Confinement has not much to do with 
it - contrary to what has been a common wisdom a decade ago and in what many 
people still believe. Nevertheless, confinement is a property of QCD, and one needs to 
understand the confinement mechanism. What can be said today is that confinement 
must be "soft": it should destroy neither the successes of the perturbative description 
of high-energy processes (no " string effects" there) nor the successes of instantons at 
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low momenta. At the moment I can think of two possible scenarios Q 

A) Instantons have a micro-structure, like merons ||; 

B) Confinement is due to monopoles which are massless because of instantons. 

Let me stress that, contrary to the case of matter where objects exist by them- 
selves, in field theory of the vacuum one has first of all to create the objects (like 
monopoles) which could bring in confinement. Therefore, only such objects can give 
a sizeable effect whose mass is effectively zero, so that one does not loose energy to 
create them. Note that in a sense instantons have zero mass since they have finite ac- 
tion, and action is mass times the (infinite) observation time. Merons like instantons 
have also finite action (if one takes care to cut them both at large and small distances), 
therefore there is nothing wrong in principle with merons. The only field-theoretical 
example of confinement we know today is the famous Polyakov's example in the 2+1 
dimensional Georgi-Glashow model |]31| , and it is an example of a meron type ||52 |. 



Moreover, a meron pair resembles one instanton ||, so a confinement mechanism 
based on merons might fit in well into the successful instanton vacuum. 

However, there is a more popular version of confinement on the market, as due 
to the Mandelstam-'t Hooft monopole condensation mechanism. According to this 
mechanism, monopole-like particles are somehow formed out of the gluon fields, they 
develop in time, interact and annihilate in pairs, and the crucial thing about them 
is that they form a quantum-mechanical condensate characterized by a macroscopic 
wave function - like the Cooper pairs of electrons in the superconductor. The con- 
finement of colour electric charges is then due to the dual Meissner effect. 

Contrary to the case of electrons which exist in abundance in any sample so that 
one needs just a small attraction to bind them into a condensate, monopoles have to 
be first of all created, and that costs energy. Therefore, this mechanism has a chance 
only if the mass of the monopoles is effectively zero. However the monopole mass can 
be estimated as M. m0 n ~ 47r/g 2 ■ {inverse size), "l/g 2 " appears here because the 
monopole should carry a unit magnetic flux. Therefore, in the weak coupling regime 
there is no chance for monopoles to condense. This is a point where instantons may 
help. 

Let us write down the effective action for gluons in the background field of one 



4 This section is based on our speculations with Victor Petrov. The reader is kindly asked not to 
be too severe to these qualitative considerations 
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instanton: 



Z 



I 



DAn exp - 




d 4 z dp dO do(p) exp 



2ir 2 p 2 



Fl{z)O ah t b + 



(13.1) 

This effective action has been first suggested in ref. || to describe the leading dipole- 
dipole interactions of instantons, then it has been re-derived in a more general form 
in ref. |34j . Later on Yung showed [ 5^j that its domain of applicability is wider than 
anticipated; recently this effective action has been used to derive the II interaction 



potential up to the next-to-next-to leading order p5 |. This effective action reproduces 
also the instanton field itself as an expension in p 2 /(x — z) 2 . For an anti-instanton 
one has to replace 77 — 77. 

For a grand canonical ensemble of Ps and 7's one writes an effective action 



rll 



E 



iV 4 



NJNJ. 



d 4 z... exp(...T]) 



N+ 



d 4 z... exp(...?7) 



N- 



DA. exp 



— ^rr — / d A z dp dO d (p) 
±9 



9„2 „2 

-2z-f-(FOfj) . - 

e 9 + e 



(13.2) 

If one expands the exponents (in the exponent) in powers of one observes 
that the linear term is zero owing to integrations over orientations, and the quadratic 
term actually corresponds to the renormalization of the gauge coupling due to the 
instanton medium: 



1 



1 



*2 




4(iV c 2 - 1)<? 4 
(2tt) 4 N 



dp d(p)(27ipY 



(13.3) 



g*(N?-l)V r J 

We see that if the packing fraction is large enough the effective coupling g* 2 for 
long range fields blows up. It means that if the instanton packing is close to some 
critical value, one does not need much energy to create a monopole in such medium; 
that is the necessary condition for their condensation. To get an accurate estimate 
of the critical density is not so easy, though. To that end one needs to have a good 
understanding of the usual perturbative renormalization of the charge by instantons: 
what is the precise argument of the running coupling constant g 2 in eq. ( |13.3j )? Using 
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the numbers obtained in ref. M we get that the density is about half that of the 
critical, but the uncertainty of this estimation is high: it could be close to the critical 
as well. 

Closing this section, I would like to mention that the monopole condensation of 
Mandelstam and 't Hooft is probably not what we in fact need - the confinement 
mechanism should be probably more subtle. The essence of that mechanism is the 
Landau-Ginzburg or the Higgs effect — but for dual Yang-Mills potentials, 't Hooft 



has elaborated it in some detail [54]]: all fields are classified in respect to the maxi- 
mal abelian U(l) x U(l) subgroup of the SU(3) colour group, and monopoles have 
magnetic charges in respect to those U(l) subgroups. If they condense all particles 
which carry electric charges in respect to those U(l) subgroups are confined. Particles 
which happen to be neutral are not, though. For example, two gluons out of eight are 
"photons" of these U(l) subgroups, so they are neutral and are not confined, instead 
they may obtain a "magnetic" mass which is of the order of the string tension, that 
is about 420 MeV, maybe up to a factor of 2 heavier. Probably such objects should 
show up as resonances in usual particle production, but we do not know of two such 
additional states. Even worse, quark- antiquark pairs belonging to the colour octet 
representation but having colour T 3 and Y zero are also neutral in respect to the U(l) 
subgroups, so they should also exist and be observable. There are two such additional 
states for each set of meson quantum numbers. There would be also five additional 
types of baryon states which are not colour singlets but which are neutral in respect 
to the both £7(1) subgroups. And of course nothing prevents monopoles themselves 
from getting into an experimentalist's detector, if only they do not, in addition, carry 
electric charges in respect to the U(l) subgroups. 

Therefore, I think that what we actually need in QCD is not condensation of 
monopoles but rather a pre-condensation, something of the kind of the Berezinsky- 
Kosterlitz-Thouless phase, characterized by large anomalous dimensions of the mono- 
pole (and probably also gluon) fields. To obtain that one also needs massless or 
effectively massless monopoles, and that is where instantons might help. 
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